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Abstract 

A search for stable and long-lived massive particles of electric charge \Q/e\ = 1 or 
fractional charges of 2/3, 4/3, and 5/3 is reported using data collected by the OPAL 
detector at LEP, at centre-of-mass energies from 130 to 209 GeV. These particles are 
assumed to be pair-produced in e + e~ collisions and not to interact strongly. No evidence 
for the production of these particles was observed. Model-independent upper limits on 
the production cross-section between 0.005 and 0.028 pb have been derived for scalar and 
spin-1/2 particles with charge ±1. Within the framework of the Constrained Minimal 
Supersymmetric Standard Model (CMSSM), this implies a lower limit of 98.0 (98.5) GeV 
on the mass of long-lived right (left)- handed scalar muons and scalar taus. Long-lived 
charged heavy leptons and charginos are excluded for masses below 102.0 GeV. For par- 
ticles with fractional charge ±2/3, ±4/3 and ±5/3, the upper limit on the production 
cross-section varies between 0.005 and 0.020 pb. All mass and cross-section limits are 
derived at the 95% confidence level and are valid for particles with lifetimes longer than 
10" 6 s. 
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1 Introduction 



Many searches for massive new particles predicted by extensions to the Standard Model (SM) 
assume that these particles decay promptly at the primary interaction vertex. Such searches 
are not sensitive to long-lived heavy particles which do not decay within the detectors. There 
exist, however, a number of models which predict such long-lived particles. For example, in the 
Constrained Minimal Supersymmetric Standard Model (CMSSM), for certain choices of the 
parameters, sleptons or charginos could be long-lived [1]. R-parity violating supersymmetric 
(SUSY) models [2] also allow for long-lived heavy particles and a fourth-generation heavy lepton 
could be stable [3]. In gauge-mediated supersymmetry, if the SUSY-breaking energy scale is 
sufficiently high [4], sleptons could be long-lived. Some models beyond the SM also predict the 
existence of particles with fractional electric charge. Previous searches for long-lived massive 
charged particles have been performed by the LEP collaborations with data taken at the Z° 
resonance [5], as well as with data taken at higher centre-of-mass energies, up to 209 GeV [6, 7]. 

This paper describes an update to a search for long-lived particles, referred to here as X ± , 
with rrix > ^z/2, and charge \Q/e\ = 1 or 2/3, pair-produced in the reaction e + e~ — > X + X~. 
This search has been described in detail in [6]. In this paper we also search for particles with 
fractional charges 4/3 and 5/3. All fractionally-charged particles are assumed to be colourless 
and non-strongly-interacting. To make the search for these particles as model independent as 
possible, only minimal calorimetric information has been used. Due to their large mass these 
particles would have anomalously high or low ionization energy loss dE/dx in the tracking 
chambers. This search is therefore primarily based on the precise dE/dx measurement provided 
by the OPAL jet chamber. The data were collected by the OPAL detector during 1995-2000, at 
centre-of-mass energies from 130 GeV to 209 GeV corresponding to a total integrated luminosity 
of 693.1 pb _1 as reported in Table 1. 

2 The OPAL Detector 

A complete description of the OPAL detector can be found in Ref. [8]. Here only a brief 
overview is given. The central detector comprised a system of tracking chambers, providing 
track reconstruction over 96% of the full solid angle 1 inside a 0.435 T uniform magnetic field 
parallel to the beam axis. It consisted of a two-layer silicon microstrip vertex detector, a high- 
precision vertex drift chamber (CV) with axial and stereo wires, a large-volume jet chamber 
and a set of z-chambers measuring the track coordinates along the beam direction. 

The jet chamber (CJ) is the most important detector for this analysis. It was divided into 
24 azimuthal sectors, each equipped with 159 sense wires. Up to 159 position and dE/dx 
measurements per track were thus possible, with a precision of a T< ^ pa 135 /im and o z pa 6 cm. 
When a track was matched with z-chamber hits and hits on the stereo wires of the vertex 
chamber (CV), the uncertainty on its z coordinate was pa 1 mm. The tracking detectors, 
located inside the magnet coil, provided a track momentum measurement with a resolution 
of a p /p pa (0.02) 2 + (0.0015 • p t ) 2 for tracks with the full number of hits (p t , in GeV, is the 
momentum transverse to the beam direction) and a resolution on the ionization energy loss 

1 The OPAL right-handed coordinate system is defined such that the z-axis is in the direction of the electron 
beam, the x-axis points towards the centre of the LEP ring, and 8 and <j> are the polar and azimuthal angles, 
defined relative to the +z- and +a;-axes, respectively. The radial coordinate is denoted by r. 
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measurement of approximately 2.8% for events with a large number of usable hits for 

dE/dx measurement [9]. 

A lead-glass electromagnetic calorimeter (ECAL) located outside the magnet coil covered 
the full azimuthal range with good hermeticity in the polar angle range of | cos#| < 0.984. The 
magnet return yoke was instrumented for hadron calorimetry covering the region | cos#| < 0.99 
and was surrounded by four layers of muon chambers. Electromagnetic calorimeters close to 
the beam axis completed the geometrical acceptance down to 24 mrad on both sides of the 
interaction point. 

The ionization energy loss dE/dx produced by a charged particle is a function of the electric 
charge Q and of fi'j = p/m, where p is the momentum and m the mass of the particle [9]. 
Figure 1 shows the distribution of dE/dx as a function of the apparent momentum, p/Q. 
Standard particles of charge ±1 (e, /i, ir, p, K) with high momentum {p > 0.1-^/s ) have 
dE/dx between 9 and 11 keV/cm. Massive particles with charge ±1 are expected to yield 
dE/dx > 11 keV/cm for high-mass values, m x > 0.36^, or dE/dx < 9 keV/cm for low-mass 
values, mx < 0.27^/s. The dE/dx measurement therefore provides a good tool for particle 
identification in these high- and low- mass regions. Massive particles with charge ±2/3 would 
have dE/dx > 11 keV/cm for high mass values, mx > 0.45^ or dE/dx < 9 keV/cm for 
low-mass values, m x < 0.35^. The expected dE/dx for massive particles of charge ±4/3 and 
±5/3 is greater than 11 keV/cm for all mass ranges. The search for massive particles with 
charge ±1/3 was not possible because the typical dE/dx deposit of these particles would be 
too close to the instrumental noise level. 

3 Monte Carlo Simulation 

Several different Monte Carlo programs were used to generate the signal process e + e~ — > X + X~. 
Signal events of the type e + e~ — > £ + £~ (i^ being a charged scalar lepton) were generated at 
several centre-of-mass energies up to 206 GeV using SUSYGEN [10]. The generated charged 
scalar leptons are not allowed to decay, therefore simulating a signal from heavy charged stable 
scalar particles. Similarly, events of the type e + e~ — > L + L~ and e + e" — > QQ, where L 1 * 1 
are stable heavy spin-1/2 leptons, and Q, Q are colourless stable heavy spin-1/2 particles with 
charge 2/3, 4/3, and 5/3, were generated at the same energies, using the EXOTIC [11] generator. 
All signal samples were generated with 1000 events per mass point with mass mx ranging from 
45 GeV to the kinematic limit for the centre-of-mass energy considered. The mass points were 
generated every 5 GeV with a finer binning of 1 GeV in the mass regions where we expect 
lower selection efficiencies. For the purpose of detector simulation and particle interactions, all 
particles were treated as heavy muons. The simulation of dE/dx in the central jet chamber of 
OPAL accounted for the charge, mass and momentum of the particle, as described above. 

The background was estimated using simulations of all Standard Model processes (two- 
fermion, four-fermion and two-photon processes) for all centre-of-mass energies from 130 to 
206 GeV. Small differences in the centre-of-mass energies of data and background Monte Carlo 
samples have a negligible effect on the analysis. 

The contribution to the background from two-fermion final states was estimated using BH- 
WIDE [12] for the e+e" final states and KORALZ [13] and KK2f [14] for the fj+fj,- and re- 
states. Hadronic two-fermion events, qq, were simulated using PYTHIA [15]. For the two- 
photon background, the PYTHIA [15], PHOJET [16] and HERWIG [17] Monte Carlo genera- 
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Figure 1: (a) The distribution of the ionization energy loss, dE/dx, of tracks in the CJ detec- 
tor, as a function of the apparent momentum, p/Q, for a sample of the data collected 
in the year 2000. The two shaded regions are the search regions. The momentum 
lower limit of the search regions is defined by the preselection cut p > 0.07^. No 
cut has been applied to the data, apart from a cutoff ofp t > 0.1 GeV made to reject 
low momentum tracks trapped in the jet chamber volume. 

(b ) Expanded view of the search regions. The theoretical curves for heavy long-lived 
particles are shown with example points from various centre-of-mass energies. In 
e + e~ — > X + X~ events, the momentum of the X particles of a given mass is fixed 
by y/s. 
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tors were used for e + e~qq final states and the Vermaseren [18] and the BDK [19] generators for 
all e + e~£ + £~ final states. Four-fermion final states were simulated with grc4f [20], which takes 
into account interference between all diagrams. All generated signal and background events 
were processed through the full simulation of the OPAL detector [21]; the same event analysis 
chain was applied to the simulated events and to the data. 

4 Data Analysis 

Pair-produced stable or long-lived massive charged particles would manifest themselves in 
events with two back-to-back tracks. Assuming they would not interact strongly, these particles 
would not produce hadronic showers. Since they are massive, they would not produce electro- 
magnetic showers either. For these reasons, the considered events would be very similar to 
events, the only difference being the higher mass of the particles, which yields a different 
dE/dx for the same momentum. 

A preselection similar to the one described in [6] is used for the analyses. Several criteria 
have been loosened in order to increase the sensitivity to high-mass particles. The criteria are 
listed below: 

PI Events are rejected if the total multiplicity of tracks in the central detector and clusters 
in the ECAL is greater than 26. Cosmic ray events are rejected as in [22]. Bhabha events 
are identified by requiring two energetic and collinear clusters in the electromagnetic 
calorimeter, these events are then rejected. 

P2 Events are required to contain at least two tracks in the central detector, each satisfying 
basic quality criteria 2 and having a momentum p > 0.07y/s, a momentum transverse to 
the beam axis p t > 0.025^, a polar angle satisfying | cos^l < 0.97 and at least 20 CJ hits 
usable for dE/dx measurement. The two selected tracks are required to have opposite 
electric charge. 

P3 To reduce background from two-photon events, the total visible energy 3 of the event 
is required to be E vis > 0.14-^/s and the acoplanarity angle 4 between the two tracks is 
required to be aC o P < 20°. 

P4 Events containing an isolated ECAL cluster with an energy greater than 5 GeV are 
rejected to reduce background from events with initial state radiation. Isolation is defined 
as an angular separation of more than 15° from the closest track. 

P5 It is required that + || < 0.2, where pi )2 are the momenta of the two selected tracks 
and E 1>2 denote the energies of the ECAL clusters associated to the tracks. This fur- 
ther reduces the contribution from Bhabha scattering events. No other tracks with 
p > 0.5 GeV and no unassociated clusters with E > 3 GeV should be found in a cone of 

2 Thc distance between the beam axis and the track at the point of closest approach (PCA) must be less 
than 1 cm; the z-coordinate of the PCA must be less than 40 cm; the innermost hit of the track measured by 
the jet chamber must be closer than 75 cm to the beam axis. 

3 The visible energy, the visible mass and the total transverse momentum of the event are calculated using 
tracks and calorimeter clusters, correcting for double counting as described in [24]. 

4 The acoplanarity angle, (f> acop , is defined as 180° minus the angle between the two tracks in the r — 4> plane. 
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10° half-opening angle around each of the two selected tracks. Criterion P5 is not used 
in the fractionally-charged particles analysis since their interaction properties with the 
calorimeters are unknown. This reduces the dependence on calorimeter response around 
the candidate tracks. 

After the preselection, the background is dominated by e + e~ — > events, with a small 

contribution from e + e~ — > t + t~ and two-photon e + e~fi + fi~ events. The effect of the preselec- 
tion cuts on the all data samples and various Monte Carlo background processes, for the search 
for \Q/e\ = 1 particles, is shown in Table 2. 



Nominal 
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1 search 


Fractional charge search 


V<GeV) 


C 


candidates 


background 


candidates 


background 


bins 


(pb- 1 ) 










133 


10.7 





0.02±0.16 





0.24±0.19 


161 


10.0 





0.11±0.11 





0.20±0.32 


172 


10.4 





0.001±0.04 





0.08±0.10 


183 


56.3 





0.13±0.33 


1 


0.37±0.95 


189 


172.3 





0.17±0.28 


1 


0.90±1.41 


192 


29.0 





0.03±0.35 





0.10±0.42 


196 


72.5 





0.14±0.40 





0.22±0.47 


200 


74.0 





0.08±0.33 





0.14±0.43 


202 


37.0 





0.00±0.18 





0.07±0.19 


205 


87.4 





0.17±0.43 





0.34±0.76 


207 


133.5 





0.26±0.66 


1 


0.52±1.16 


Total 


693.1 





1.1±1.3 


3 


3.2±2.4 



Table 1: The number of candidate events and the expected background at all energies, for 
the search for \Q/e\ = 1 and fractionally-charged particles. The errors quoted in- 
clude both statistical and systematic effects. In the second column, the integrated 
luminosity is given for each energy. The data collected in the year 2000 were de- 
livered at various centre-of-mass energies, up to yfs = 209 GeV. For this analysis 
they have been separated in two bins, the first one, referred to as 205 GeV, with 
y/s < 206 GeV and an average yfs of 204. 7 GeV, and the second one, referred to as 
207 GeV, with yfi > 206 GeV and an average yfs of 206.6 GeV. 



4.1 Search for particles with unit charge 

The search strategy has been simplified with respect to [6] and now relies entirely on dE/dx 
information. A sample dE/dx distribution for data and simulated Monte Carlo events is shown 
in Figure 2. The kinematic selection is no longer used. This new strategy has been applied to 
all data samples. The preselected events are retained if they satisfy the following requirements 
on dE/dx: 

Al Both high- momentum tracks must have either dE/dx > 11 keV/cm or dE/dx < 9 
keV/cm. 
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Figure 2: dE/dx distribution for data and Monte Carlo simulation at y/s = 189 GeV after 
the preselection cut P4. The arrows show the accepted region. 



A2 The probability that the dE/dx measurements for either track were consistent with one 
of the standard particles (e, //, ir, p, K) must be less than 30%. This removes background 
from poorly measured SM particles. 



Cuts 


Data 


Background Simulation 


Signal MC (%) 


Total 


e + e 






e + e~£ + £- 


Others 






e 80 


Pl-2 


20199 


19572.0 


9292.8 


4015.2 


1305.4 


3428.9 


1529.7 


96.6 


97.5 
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3327.2 
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93.0 


94.3 


P5 
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4956.6 


<1.0 


3098.7 


76.5 


1671.3 


110.1 
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1.1 


<1.0 


0.01 


0.01 


1.1 


0.03 


85.7 


73.4 


94.3 



Table 2: For the \Q/e\ = 1 analysis, the numbers of events remaining after each cut for all 
data collected at the various centre- of- mass energies and for various Monte Carlo 
background processes normalised to the integrated luminosity of the data ("Others" 
refers to e + e~ — > qq and e + e~ — > four-fermion processes) . When no candidate events 
are selected in the Monte Carlo, a 68% CL upper limit on the number of events is 
used as the statistical uncertainty. In the last three columns, the efficiencies for 
are given (in percent) for mx = 45, 55, 80 GeV at *Js = 207 GeV. 

The effect of cuts A1-A2 for all data and simulated events can be seen in Table 2. No 
candidate event is found. The expected backgrounds at the various centre-of-mass energies are 
shown in Table 1. The total background, summed over all energies, is estimated to be 1.1±1.3 
events, where the uncertainty quoted includes both statistical and systematic effects. 

The detection efficiency for spin-0 particles varies between 75 and 90% for masses mx < 
0.27\/s or mx > 0.36^. The efficiency drops significantly in the region 0. 27^/5 < mx < 
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0.36^, but as data sets collected at different centre-of-mass energies (y/s) are combined, a 
reasonable selection efficiency is achieved for all mass values up to close to the kinematic limit. 
The selection efficiencies of all Monte Carlo samples at all y/s are parametrised as a function 

of /?7 = p/mx = \J~sj Am? x — 1 of the particle. This parametrisation is used to calculate the 
efficiency for all masses at each centre-of-mass energy, using linear interpolation. For spin-1/2 
particles, the efficiencies are 2-9% lower than for spin-0 particles due to the different angular 
distribution of the tracks. We analyse each centre-of-mass energy separately, then combine the 
results for the final cross-section limits. 

4.2 Search for particles with fractional charges 

To search for particles with fractional charges of 2/3, 4/3, and 5/3 the selection criteria PI 
through P4, followed by Al and A2 are used. The results after each cut for all data samples and 
various Monte Carlo background processes are shown in Table 3. For charge 2/3 the selection 
efficiency is between 75 and 90% for most of the mass range, while for charges 4/3 and 5/3 the 
efficiency is above 90% for the whole mass range. After this selection, one candidate survives in 
the data sample at yfs = 183 GeV, one at yfs = 189 GeV, and one at y/s = 207 GeV, while no 
candidate is left in any of the other data sets. The masses of the candidates are reconstructed 
for charge 2/3 using the dE/dx and momentum information, while for charge 4/3 and 5/3 
kinematic information only is used. The reconstructed masses, track momenta and dE/dx 
values of the candidate events are reported in Table 4. The selected events and the expected 
background at each centre-of-mass energy are shown in Table 1. The total background, summed 
over all energies, is estimated to be 3.2±2.4 events, where the uncertainty quoted includes both 
statistical and systematic effects. 



Cuts 


Data 


Background Simulation 


Signal MC (%) 






Total 


e + e 






e + e~£ + £- 


Others 


2/3 


4/3 


5/3 


Pl-2 


20199 


19572.0 


9292.8 


4015.2 


1305.4 


3428.9 


1529.7 


93.1 


97.6 


97.8 


P3-4 


15935 


15447.4 


7439.8 


3327.2 


1081.0 


2952.8 


646.6 


89.7 


94.8 


94.7 


A1,A2 


3 


3.2 


0.7 


0.03 


0.5 


1.4 


0.5 


86.3 


94.8 


94.7 



Table 3: For fractional charges analyses, the numbers of events remaining after each cut for 
all data collected at the various centre-of-mass energies and for various Monte Carlo 
background processes normalised to the integrated luminosity of the data ("Others" 
refers to e + e~ — > qq and e + e~ — > four-fermion processes) . In the last three columns, 
the efficiencies are given (in percent) for \Q/e\ = 2/3, 4/3, 5/3 at mx = 70 GeV and 
v/i = 207 GeV. 



4.3 Systematic uncertainties 

The main sources of systematic uncertainties for this analysis are discussed below and reported 
in Table 5, for both the signal efficiency and the background estimate: 

• The errors on the Monte Carlo modeling of acO p, -^vis and dE/dx are estimated by 
comparing the distributions of these variables for data and background Monte Carlo. 
The relative difference between the averages of the distributions is used to increase or 
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Pi 


Pi 


(dE/dx) i 


(dE/dx) 2 


Masses (GeV) 


(GeV) 


(GeV) 


(GeV) 


(keV/cm) 


(keV / cm) 


\QM = 1 


\QM = § 


\QM = 1 


183 


64 ± 15 


19 ±30 


11.04 ±0.62 


20.48 ± 1.00 


49.6 ± 11.9 


95.9 ±9.4 




189 


30 ±2 


29 ± 15 


11.17 ±0.39 


19.85 ±0.95 


24.9 ± 1.3 


85.9 ± 1.0 


80.7 ± 1.6 


207 


77 ±8 


71 ±8 


11.06 ±0.34 


11.13 ±0.36 


61.7 ±4.7 







Table 4: Information on the candidate events selected by the fractionally-charged analysis. 

The momentum and dE/dx of each track is reported together with the reconstructed 
masses for the \Q/e\ = 2/3, 4/3, 5/3 hypothesis. The candidate mass is not reported 
when the reconstruction procedure gives a kinematically inconsistent result. 



Quantity 


Systematic uncertainty (%) 


Signal 


Background 


High efficiency region 


Low efficiency region 


0acop 

E ■ 

■'-'VIS 

dE/dx 
MC statistics 
Interpolation 
MC generator 
Double tracks 
Luminosity 


0.0-1.2 
0.0-2.1 


0.0 

0.0 
<0.008 
58-157 


0.0-0.3 
0.6-2.6 
0.0-2.7 


0.0-27 
0.8-21 
0.0-35 




5.9-300 
2.9 


0.22-0.68 


Total 


0.7-3.6 


1.3-39 


90-309 



Table 5: Relative systematic uncertainties in the signal efficiency associated with the various 
quantities used for the \Q/e\ = 1 search. The ranges given cover all centre-of-mass 
energies. The systematic errors vary slightly with centre-of-mass energy, but strongly 
with mx for the signal. For this reason the two regions are reported: high efficiency 
(mx/y/s < 0.27 or mx/y/s > 0.36) and low efficiency (0.27 < mx/y/s < 0.36). 

decrease the value of the cut on the relevant variable in order to decrease the overall 
signal efficiency and background estimate. The difference between the reduced efficiency 
and the one obtained with the nominal selection is taken as the systematic uncertainty 
due to the modeling of the variable under consideration. This estimate of the error is 
more conservative than the one obtained by smearing the dE/dx values of the tracks. 

• The MC statistical uncertainty, due to the limited number of signal events generated, has 
been computed using a binomial formula. For background processes the large relative 
statistical uncertainty is due to the limited number of background events selected. 

• The uncertainty due to the linear interpolation of the signal detection efficiency is esti- 
mated as the difference between the interpolated values and the efficiency obtained at 
mass points where MC signal samples were generated, when that mass point was omitted 
from the interpolation procedure. 

• For background processes of the type e + e~£ + £~ , the Vermaseren [18] generator has been 
used as the reference generator. The difference in the background expectation obtained 
by using BDK [19] instead of Vermaseren has been considered as the uncertainty on the 
MC generator. 
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• Two tracks separated by a distance smaller than 2.5 mm could be unresolved and recon- 
structed as a single track with a high dE/dx value. Events with unresolved double tracks 
are potential backgrounds for this search. The systematic uncertainty introduced by the 
modeling of the double track resolution in the Monte Carlo samples has been estimated 
to be 2.9 percent. 

The absolute uncertainty on the background is reported by centre-of-mass energy in Ta- 
ble 1. The uncertainty introduced on the integrated luminosity [25] is also reported in Table 5. 
At a given centre-of-mass energy the different systematic uncertainties are assumed to be in- 
dependent, so that the total systematic uncertainty is calculated as the quadratic sum of the 
individual uncertainties. The Monte Carlo modeling of dE/dx, the Vermaseren-BDK genera- 
tor, and luminosity uncertainties are assumed to be correlated between centre-of-mass energies, 
while the other systematic uncertainties are assumed to be independent. 

5 Results 

The numbers of candidates found in the search for particles with charge \Q/e\—l and fractional 
charges are summarised for all energies in Table 1, together with the expected backgrounds. 
The data show no significant excess above the expected background from Standard Model 
processes. 

Model-independent cross-section upper limits have been computed for the pair-production 
of massive charged long-lived particles, combining the results from all centre-of-mass energies, 
assuming s-channel production. The cross-section dependence on the energy is taken to be f3 z / s 
for spin-0 particles and |(1 — ! y) for spin-1/2 particles, where (5 = p/E ~ \A — 4m|/s. In 
evaluating upper limits, the candidates are counted in mass intervals centred on their central 
mass values and ±2<r wide (where a is the error on the mass estimate of the candidates as 
reported in Table 4). In the case in which the mass could not be reconstructed the candidates 
were considered in the whole mass range (from 45 GeV to the kinematic limit for that event). 
A likelihood-ratio method [26] was used to determine an upper limit for the cross-section. The 
total systematic error is incorporated into the limits following the prescription of Ref. [27]. 

In Figure 3(a), the 95% C.L. upper limit on the cross-section at \fs = 206.6 GeV is shown 
for spin-0 particles of charge ± 1. The 95% C.L. upper limit on the pair-production cross- 
section varies from 0.005 to 0.028 pb in the mass range 45 < mx < 101 GeV. The bump in 
the mass range of 52 < m x < 70 GeV is due to the low efficiencies described in section 4.1. 
The cross-section limits are compared with the predicted cross-sections [10] for pair-production 
of right- and left-handed smuons and staus to determine mass limits. For these two slepton 
species, the production cross-section does not depend on the CMSSM parameters but only on 
the slepton mass. The 95% C.L. lower mass limits are 98.0 GeV and 98.5 GeV for the mass of 
right- and left-handed smuons and staus, respectively, as shown in Figure 3(a). 

Figure 3(b) shows the 95% C.L. upper limit on the cross-section at y/s = 206.6 GeV for 
spin-1/2 particles of charge ±1. The limit varies from 0.005 to 0.024 pb in the mass range 
45 < m x < 100 GeV. This limit must be compared with the predicted cross-sections for 
chargino production [10] and heavy charged lepton production [11, 28]. For the chargino limits, 
the CMSSM parameters have been chosen to minimise the predicted chargino cross-section 
at every chargino mass value (assuming a heavy sneutrino, rriy > 500 GeV), without any 
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Figure 3: Model-independent 95% C.L. upper limits on the pair-production cross-section of 
spin-0 (a) and spin- 1/2 (b) heavy long-lived non-strongly-interacting particles of 
charge ±2 as a function of their mass (solid line) at \fs = 206.6 GeV. The bump 
observed between masses of 52 and 70 GeV is due to the drop in efficiency where the 
dE/dx expected for signal crosses the band of dE/dx values expected for standard 
particles. The CMSSM predicted cross-sections for right-handed (dash-dotted line) 
and left-handed (dashed line) smuons and staus are also shown in (a). The 95% 
C.L. lower limits on the masses of these sleptons are at the crossing point between 
the experimental limit and theoretical prediction. 

restriction on the mass of the lightest neutralino. A 95% C.L. lower limit on the masses of 
long-lived charginos, of 102.0 GeV, is obtained for every choice of the CMSSM parameters. 
The 95% C.L. lower limit on the heavy charged lepton mass is also 102.0 GeV. 

Figure 4(a) shows the 95% C.L. upper limit on the cross-section at y/s = 206.6 GeV for spin- 
1/2 particles of charge ±2/3. The limit varies between 0.005 and 0.020 pb in the mass range 
45 < mx < 101 GeV. Figures 4(b) and 4(c) show the 95% C.L. upper limit on the cross-section 
at a/s = 206.6 GeV for spin-1/2 particles of charge ±4/3 and of charge ±5/3, respectively. For 
spin-0 particles with fractional charge the cross-section upper limits are slightly higher than for 
spin-1/2 particles. This is due to the difference in angular distributions. 

All results obtained are valid for non-strongly-interacting colourless particles with a lifetime 
longer than 10~ 6 s. This lifetime restriction is obtained by considering the heaviest (and there- 
fore slowest) particles excluded by this search, and then requiring that the decay probability 
of these particles at a flight distance larger than 3.0 m be greater than 95%. For lower mass 
values the results are also valid for shorter lifetimes. 



6 Summary and Conclusions 

A search was performed for pair-production of stable and long-lived massive particles not sub- 
ject to strong interactions, with charge \Q/e\ = 1 or fractional charges of 2/3, 4/3, and 5/3. The 
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Figure 4: Model-independent 95% C.L. upper limits on the pair-production cross-section as 
a function of mass of spin- 1/2 heavy long-lived non-strongly-interacting particles 
of charge (a) \Q/e\ = ±2/3, (b) \Q/e\ = ±4/3 and (c) \Q/e\ = ±5/3, at ^1 = 
206.6 GeV. 

primary tool used in this search was the precise dE/dx measurement provided by the OPAL 
jet chamber. No evidence for the production of such particles was observed. For s-channel 
production, the upper limits on the cross-section vary between 0.005 and 0.026 pb in the mass 
range explored for particles of charge ±1. Within the framework of the CMSSM, lower mass 
limits have been obtained: 98.0 GeV for right-handed and 98.5 GeV for the left-handed smuons 
and staus. Charged long-lived massive leptons and long-lived charginos with masses smaller 
than 102.0 GeV are excluded. For particles with fractional charge ±2/3, ±4/3 and ±5/3, the 
upper limits on the production cross-section vary between 0.005 and 0.020 pb in the range 
45 < mx < 95 GeV. The above limits are valid at the 95% C.L. for particles with lifetimes 
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longer than 10 6 s. 



Acknowledgements 

We particularly wish to thank the SL Division for the efficient operation of the LEP accelerator 

at all energies and for their close cooperation with our experimental group. In addition to the 

support staff at our own institutions we are pleased to acknowledge the 

Department of Energy, USA, 

National Science Foundation, USA, 

Particle Physics and Astronomy Research Council, UK, 

Natural Sciences and Engineering Research Council, Canada, 

Israel Science Foundation, administered by the Israel Academy of Science and Humanities, 
Benoziyo Center for High Energy Physics, 

Japanese Ministry of Education, Culture, Sports, Science and Technology (MEXT) and a grant 

under the MEXT International Science Research Program, 

Japanese Society for the Promotion of Science (JSPS), 

German Israeli Bi-national Science Foundation (GIF), 

Bundesministerium fur Bildung und Forschung, Germany, 

National Research Council of Canada, 

Hungarian Foundation for Scientific Research, OTKA T-038240, and T-042864, 
The NWO/NATO Fund for Scientific Reasearch, the Netherlands. 



References 

[1] H. P. Nilles, Phys. Rep. 110 (1984) 1; 

H. E. Haber and G. L. Kane, Phys. Rep. 117 (1985) 75. 

[2] H. Dreiner, "An Introduction to Explicit R-parity Violation", Perspective on Supersym- 
metry, ed. G.L Kane, 462, and hep-ph/9707435. 

[3] E. Nardi, E. Roulet and D. Tommasini, Phys. Lett. B344 (1995) 225. 

[4] S. Dimopoulos, S. Thomas and J. D. Wells, Nucl. Phys. B488 (1997) 39; 
S. Ambrosanio, G. D. Kribs and S. P. Martin, Phys. Rev. D56 (1997) 1761; 
G. F. Giudice and R. Rattazzi, Phys. Rep. 322 (1999) 419. 

[5] ALEPH Collab., D. Decamp et al, Phys. Lett. B303 (1993) 198; 
DELPHI Collab., P. Abreu et al, Phys. Lett. B247 (1990) 157; 
OPAL Collab., K. Ahmet et al, Phys. Lett. B252 (1990) 290; 
OPAL Collab., R. Akers et al, Z. Phys. C67 (1995) 203. 

[6] OPAL Collab., R. Ackerstaff et al, Phys. Lett. B433 (1998) 195. 

[7] ALEPH Collab., R. Barate et al, Phys. Lett. B405 (1997) 379; 
DELPHI Collab., P. Abreu et al, Phys. Lett. B478 (2000) 65; 



15 



DELPHI Collab., P. Abreu et al, Phys. Lett. B503 (2001) 34; 
L3 Collab., P. Achard et al, Phys. Lett. B517 (2001) 75. 

[8] OPAL Collab., K. Ahmet et al, Nucl. Instr. Meth. A305 (1991) 275; 
P. P. Allport et al, Nucl. Instr. Meth. A324 (1993) 34; 
P. P. Allport et al, Nucl. Instr. Meth. A346 (1994) 476; 
S. Anderson et al, Nucl. Instr. Meth. A403 (1998) 326. 

[9] M. Hauschild et al, Nucl. Instr. Meth. A314 (1992) 74; 
M. Hauschild, Nucl. Instr. Meth. A379 (1996) 436. 

[10] S. Katsanevas and S. Melachroinos, in "Physics at LEP2", ed. G. Altarelli, T. Sjostrand 
and F. Zwirner, CERN 96-01, vol. 2, (1996) 328; 

S. Katsanevas and P. Morzwitz, Comp. Phys. Comm. 112 (1998) 227. 
[11] R. Tafirout and G. Azuelos, Comp. Phys. Comm. 126 (2000) 244. 

[12] S. Jadach et al, in "Physics at LEP2", eds. G. Altarelli, T. Sjostrand and F. Zwirner, 
CERN 96-01, vol.2, (1996) 229; 

S. Jadach, W. Placzek and B.F.L. Ward, Phys. Lett. B390 (1997) 298. 

[13] S. Jadach, B.F.L. Ward and Z. Was, Comp. Phys. Comm. 79 (1994) 503. 

[14] S. Jadach, B.F.L. Ward and Z. Was, Comp. Phys. Comm. 130 (2000) 260. 

[15] T. Sjostrand and M. Bengtsson, Comp. Phys. Comm. 43 (1987) 367; 

"PYTHIA 5.7 and JETSET 7.4, Physics and Manual", CERN-TH. 7112/93 (revised Au- 
gust 1995); 

T. SjSstrand, Comp. Phys. Comm. 82 (1994) 74; 
T. Sjostrand, Comp. Phys. Comm. 135 (2001) 238. 

[16] R. Engel and J. Ranft, Phys. Rev. D54 (1996) 4244; 
R. Engel, Z. Phys. C66 (1995) 203. 

[17] G. Marchesini et al, Comp. Phys. Comm. 67 (1992) 465. 

[18] R. Bhattacharya, J. Smith and G. Grammer, Phys. Rev. D15 (1977) 3267; 
J. Smith, J. A.M. Vermaseren and G. Grammer, Phys. Rev. D15 (1977) 3280; 
J. Smith, J. A.M. Vermaseren and G. Grammer, Phys. Rev. D19 (1979) 137. 

[19] F.A.Berends, PH. Daverveldt and R. Kleiss, Nucl. Phys. B253 (1985) 421; 

F.A.Berends, P.H. Daverveldt and R. Kleiss, Comp. Phys. Comm. 40 (1986) 271; 
F.A.Berends, P.H. Daverveldt and R. Kleiss, Comp. Phys. Comm. 40 (1986) 285; 
F.A.Berends, P.H. Daverveldt and R. Kleiss, Comp. Phys. Comm. 40 (1986) 309. 

[20] J. Fujimoto et al, Comp. Phys. Comm. 100 (1997) 128. 

[21] J. Allison et al, Nucl. Instr. Meth. A317 (1992) 47. 

[22] OPAL Collab., R. Akers et al, Z. Phys. C61 (1994) 19. 

[23] S. Robins, "A Study of Bhabha Scattering at LEP", Ph.D. Thesis (1991), Queen Mary, U. 
of London, RALT-136. 



16 



[24] OPAL Collab., M.Z. Akrawy et al, Phys. Lett. B253 (1991) 511. 
[25] OPAL Collab., K. Ackerstaff et al, Phys. Lett. B391 (1997) 221. 

[26] A. G. Frodesen, O. Skeggestad and T. Tofte, "Probability and Statistics in Particle 
Physics", Universitetsforlaget, 1979, ISBN 82-00-01-01906-3; 

S. L. Meyer, "Data Analysis for Scientists and Engineers", John Wiley and Sons, 1975, 
ISBN 0-471-59995-6. 

[27] R. D. Cousins and V. L. Highland, Nucl. Instr. Meth. A320 (1992) 331. 

[28] A. Djouadi, Z. Phys. C63 (1994) 317. 



17 



